Spectral Characteristic Evolution: A New Algorithm for Gravitational Wave Propagation

We present a spectral algorithm for solving the full nonlinear vacuum Einstein field equations in the Bondi framework. Developed within the Spectral Einstein Code (SpEC), we demonstrate spectral characteristic evolution as a technical precursor to Cauchy Characteristic Extraction (CCE), a rigorous method for obtaining gauge-invariant gravitational waveforms from existing and future astrophysical simulations. We demonstrate the new algorithm's stability, convergence, and agreement with existing evolution methods. We explain how an innovative spectral approach enables a two orders of magnitude improvement in computational efficiency.Comment: 28 pages, 9 figure

Key elements of robustness in binary black hole evolutions using spectral methods

As a network of advanced-era gravitational wave detectors is nearing its design sensitivity, efficient and accurate waveform modeling becomes more and more relevant. Understanding of the nature of the signal being sought can have an order unity effect on the event rates seen in these instruments. The paper provides a description of key elements of the Spectral Einstein Code (SpEC), with details of our spectral adaptive mesh refinement (AMR) algorithm that has been optimized for binary black hole (BBH) evolutions. We expect that the gravitational waveform catalog produced by our code will have a central importance in both the detection and parameter estimation of gravitational waves in these instruments

Spectral Cauchy Characteristic Extraction of strain, news and gravitational radiation flux

We present a new approach for the Cauchy-characteristic extraction of gravitational radiation strain, news function, and the flux of the energy-momentum, supermomentum and angular momentum associated with the Bondi-Metzner-Sachs asymptotic symmetries. In Cauchy-characteristic extraction, a characteristic evolution code takes numerical data on an inner worldtube supplied by a Cauchy evolution code, and propagates it outwards to obtain the space-time metric in a neighborhood of null infinity. The metric is first determined in a scrambled form in terms of coordinates determined by the Cauchy formalism. In prior treatments, the waveform is first extracted from this metric and then transformed into an asymptotic inertial coordinate system. This procedure provides the physically proper description of the waveform and the radiated energy but it does not generalize to determine the flux of angular momentum or supermomentum. Here we formulate and implement a new approach which transforms the full metric into an asymptotic inertial frame and provides a uniform treatment of all the radiation fluxes associated with the asymptotic symmetries. Computations are performed and calibrated using the Spectral Einstein Code (SpEC).Comment: 30 pages, 17 figure

Batch Cooling Crystallization of Plate-like Crystals: A Simulation Study

In the paper, batch cooling crystallization of plate-like crystals is investigated by numerical experimentation. The population balance approach is used to describe a seeded crystallization process, developing a bivariate morphological population balance model including nucleation and size dependent crystal growth. The population balance equation is reduced to a system of equations for the bivariate mixed moments of internal variables, which are solved by means of the quadrature method of moments. A sensitivity analyses is carried out to investigate the effects of the process parameters on dynamic evolution of particulate properties including the characteristic mean particle size and shape

Joint approach for reducing eccentricity and spurious gravitational radiation in binary black hole initial data construction

At the beginning of binary black hole simulations, there is a pulse of spurious radiation (or junk radiation) resulting from the initial data not matching astrophysical quasi-equilibrium inspiral exactly. One traditionally waits for the junk radiation to exit the computational domain before taking physical readings, at the expense of throwing away a segment of the evolution, and with the hope that junk radiation exits cleanly. We argue that this hope does not necessarily pan out, as junk radiation could excite long-lived constraint violation. Another complication with the initial data is that they contain orbital eccentricity that needs to be removed, usually by evolving the early part of the inspiral multiple times with gradually improved input parameters. We show that this procedure is also adversely impacted by junk radiation. In this paper, we do not attempt to eliminate junk radiation directly, but instead tackle the much simpler problem of ameliorating its long-lasting effects. We report on the success of a method that achieves this goal by combining the removal of junk radiation and eccentricity into a single procedure. Namely, we periodically stop a low resolution simulation; take the numerically evolved metric data and overlay it with eccentricity adjustments; run it through an initial data solver (i.e. the solver receives as free data the numerical output of the previous iteration); restart the simulation; repeat until eccentricity becomes sufficiently low; and then launch the high resolution “production run” simulation. This approach has the following benefits: (1) We do not have to contend with the influence of junk radiation on eccentricity measurements for later iterations of the eccentricity reduction procedure. (2) We reenforce constraints every time the initial data solver is invoked, removing the constraint violation excited by junk radiation previously. (3) The wasted simulation segment associated with the junk radiation’s evolution is absorbed into the eccentricity reduction iterations. Furthermore, (1) and (2) together allow us to carry out our joint-elimination procedure at low resolution, even when the subsequent “production run” is intended as a high resolution simulation

Improvements to the construction of binary black hole initial data

Construction of binary black hole initial data is a prerequisite for numerical evolutions of binary black holes. This paper reports improvements to the binary black hole initial data solver in the Spectral Einstein Code, to allow robust construction of initial data for mass-ratio above 10:1, and for dimensionless black hole spins above 0.9, while improving efficiency for lower mass-ratios and spins. We implement a more flexible domain decomposition, adaptive mesh refinement and an updated method for choosing free parameters. We also introduce a new method to control and eliminate residual linear momentum in initial data for precessing systems, and demonstrate that it eliminates gravitational mode mixing during the evolution. Finally, the new code is applied to construct initial data for hyperbolic scattering and for binaries with very small separation.Comment: 28 pages, 13 figures, 1 tabl